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Abstract -The analysis of stability of community structure is an important problem for scientists 
from many fields. Here, we propose a new framework to reveal hidden properties of community 
structure by quantitatively analyzing the dynamics of Potts model. Specifically we model the Potts 
procedure of community structure detection by a Markov process, which has a clear mathematical 
explanation. Critical topological information regarding to multivariate spin configuration could 
also be inferred from the spectral significance of the Markov process. We test our framework on 
some example networks and find it doesn’t have resolute limitation problem at all. Results have 
shown the model we proposed is able to uncover hierarchical structure in different scales effectively 
and efficiently. 
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j 1. Introduction. — Community structure detection 
O00 is a main focus of complex network studies. It has 
00 attracted a great deal of attentions from various scientific 
<0 fields. Intuitively, community refers to a group of nodes 
CO in the network that are more densely connected internally 
<0>than with the rest of the network. A well known explo- 
ration for this problem is the concept of modularity, which 
Tl is proposed by Newman et al mm to quantify a network’s 
L* partition. Optimizing modularity is effective for commu- 
nity structure detection and has been widely used in many 
real networks. However, as pointed out by Fortunato et al 
cm modularity suffers from the resolution limit problem 
which is concerned about the reliability of the communities 
detected through the optimization of modularity. Com¬ 
plementary to the modularity concept, many efforts are 
devoted to understanding the properties of the dynamical 
processes taking place in the underlying networks. Specif¬ 
ically, researchers have begun to investigate the correla¬ 
tion between the community structure and the dynamical 
systems, such as synchronization s 0 and random walk 
process [B] [15]. 
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Potts dynamical model is a powerful tool which has been 
applied to uncover the thermodynamical behaviors in net¬ 
works 00. It models an inhomogeneous ferromagnetic 
system where each node is viewed as a labeled spin in the 
network. The configuration of the system is defined by 
the interactions between the nodes. Considering an un¬ 
weighted network with N nodes without self-loops, a spin 
configuration {5} is defined by assigning each node i a spin 
label Si which may take integer values Si = 1, • • ■, K. To 
characterize the coherence between two nodes, spin-spin 
correlation C = (Cij) is defined as the thermal average of 

S Si , Sj - 

= (S Si , Sj ), (1) 

which represents the probability that spin variables s„; and 
Sj have the same value. Cij takes values from the inter¬ 
val [0,1], representing the continuum from no coupling to 
perfect accordance of nodes i and j. In section 4, we de¬ 
velop a novel hierarchical block model which can calculate 
C efficiently. The C is corresponding to the average spin 
correlation across multiple level of the hierarchical struc¬ 
ture. 

If the system is not homogeneous but has a community 
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time(r) 



Fig. 1: Dynamics of spin configuration of four communities 
when they go through several local uniform states to the global 
stable state. Different spin values are described by different 
shapes in dynamical tree, n represents the time at which the 
system has i different spin states 

. The significance curve of spin configurations at 
different times can be calculated shown in the left part of 
the graph. 


network. Markov process is a useful tool and has been 
applied to find communities [B]. Let P = ( pij ) be the 
stochastic transition matrix and the element p^ is defined 
as 


Ct. 


Pij — 


Ef=i Cij 


( 2 ) 


where Cij is the spin-spin correlation function defined in 
Eq.{IJ). Via this representation, the tools of stochastic 
theory and finite-state Markov processes [ 6 ] can be utilized 
for the purpose of community structure analysis. 

For this ergodic Markov process, P* represents the tran¬ 
sition probability matrix between nodes over a period 
of t time steps. To compute the transition matrix P 4 , 
the eigenvalue decomposition of P is used. If Xk with 
k = 1, • ■ •, n denote the eigenvalues of P, and its right and 
left eigenvectors Uk and fk are scaled to satisfy 


structure, the states are not just ferromagnetic or para¬ 
magnetic [7]. We assume that the spins will go through 
a hierarchy of local uniform states (meta-stable states) as 
time increases which is shown in Fig[lJ before they reach 
a globally stable state with the same value. In each lo¬ 
cal uniform state, spin values of nodes within the same 
communities are identical. Correspondingly, one can cal¬ 
culate the hitting and exiting time of each local uniform 
state and there should be a big gap between them when a 
well-formed community structure exists. The significance 
of spin configurations at different time can also be calcu¬ 
lated to illustrate the amplitude of variation. 

In this letter, using the Potts model and spectral the¬ 
ory, we firstly uncover the relationship between commu¬ 
nity structure of a network and its meta-stability of spin 
dynamics, and then propose the significance of communi¬ 
ties to characterize and analyze the underlying spin con¬ 
figuration. For any given network, one can straightfor¬ 
wardly get critical information related to its community 
structure, such as the stability and the optimal number 
of communities across multiple timescales without using 
particular partition algorithms. We then use phase tran¬ 
sition of stochastic dynamical system to prove that the 
stability we proposed is able to indicate the significance of 
community structure more theoretically which is based on 
eigengap theory. Furthermore, a novel hierarchical block 
model is proposed which can calculate spin correlation at 
each layer of the network structure. Finally, we test our 
framework on some examples of complex networks. Re¬ 
sults show the model we proposed is able to uncover the 
hierarchical structure in different scales effectively and ef¬ 
ficiently and doesn’t have resolute limitation problem at 
all. 

2. The framework. — In order to establish the con¬ 
nection between the community structure and the local 
uniform behavior of Potts model, we introduce a Markov 
stochastic model featured by spectral significance for the 


Puk = XkUk , fkP = Xkfk ( 3 ) 

The orthonormality relation of Uk and fi is satisfied: 

u k fi = Ski, (4) 

and the spectral representation of P is given by 

i* = (5) 

k 

We assume that eigenvalues of P are sorted such that 
Ai = 1 > |A 2 1 > |A 3 1 > ... > |A„|. Because of the largest 
eigenvalue Ai = 1, when time t —» 00 , P = P°° = uifi. 
The convergence of every initial distribution to the sta¬ 
tionary distribution P^ corresponds to the fact that the 
spins of whole system ultimately reach exactly the same 
value, as time increases. This perspective belongs to a 
timescale t —» 00 , at which all eigenvalues A), go to 0 ex¬ 
cept for the largest one, A) = 1. In the other extreme 
of a timescale t = 0, P 4 becomes the stationary distri¬ 
bution matrix. All of its columns are different, and the 
system disintegrates into as many spin values as possi¬ 
ble. Then, we simply extend P 4 to the symmetrical form 
G (i) = (pt + (pi)T)/2. 

Suppose the partition method divides the network A 
into K communities or sets 14 C V,k € {1,2 
which are disjoint and the sets 14, 14,..., Vk together 
form a partition of node set V. The number of nodes 
in each community is denoted by A4 = |I4|- We take the 
time series into consideration. Therefore, we define the 
significance of a given community k by the ratio of inner 
correlations as 


S 


(t) 

k 


E 

*>ieVfc 


[G (t) kj 

N k 


( 6 ) 


can be viewed as a function of timescale t and we can 
use it to study the trend of community structure as time 
goes on. 
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Further discussion is facilitated by reformulating the av¬ 
erage association objective in matrix form. We denote the 
membership vector of community k by x k , a binary vector 
that describes each node’s involvement in community k. 
The hard partition and disjointness of sets 14 requires that 
the vectors xi and x s are orthogonal. Given the number of 
communities K, the communities are found by maximizing 
the objective function 


= E e 


fc=l i,jGVk 


N k 


vk xlG^Xk 

h x k x fc 


(7) 


The objective is to be maximized under the conditions 
Xk £ {0,1} and xjx s = 0 if l ^ s. Eq.([7J can be rewritten 
as a matrix trace by accumulating the vectors u k into a 
matrix X = (xi,X 2 , ■■■, xx)- We can then write the objec¬ 
tive as 


= tr{(X T X)- 1 X T G^X} , . 

= tr{{X T X)- 1 / 2 X T G^X(X T X)- 1 / 2 } [ ’ 

where matrix X T X is diagonal. The substitution Y = 
X(X T X)- 1 / 2 simplifies the optimization problem to 
Jpp = tr{Y T G^Y}. The condition Y T Y = lx is au¬ 
tomatically satisfied since 


Y t Y = (X T X)- 1/2 (X T X)(X T X)~ 1/2 = I K . (9) 

The vectors y k thus have unit length and are orthogonal 
to each other. The optimization problem can be written 
in terms of the matrix Y as 


From another perspective, because the eigenvalues are 
sorted by Ai = 1 > |A 2 1 > |A 3 1 > ... > |A„|, the strength 
of a community at time t, X k , can also be viewed as the 
robustness of fc-spin state at time t. At this point, the 
eigengap X t k _ 1 — X k can be interpreted as the “difficulty” 
that the fc-spin state transfer to the (k — l)-spin state at 
time t. The number of communities A at time t is then 
inferred from the location of the maximal eigengap, and 
this maximal value can be used as a quality measure for 
the most stable state. The A (t) is formally defined as 

A(f) = arg[maxk(X t k _ 1 - A^)] (12) 

From a global perspective if the number of communities A 
doesn’t change for the longest time, we can consider it as 
the optimal number for this network, represented as ’F. 

To a certain extent, the most stable state can represent 
the spin configuration of the whole network. Thus, we de¬ 
fine the stability of community structure at each timescale, 
0 (i), as the stability of the most stable spin state: 

©(*) = ^A(t )—1 ~ ^A(t) (13) 

Our expectation is that from the trend of O (f), one can 
find the most stable timescale for community structure 
where O(i) reaches the maximal. Furthermore, from a 
global perspective, we can use the largest stability corre¬ 
sponding to q communities, T(q) = max{Q(t)\A(t) = q}, 
to indicate the robustness of a network, defined as the 
stability of the structure with q communities. While T(q) 
tries to directly characterize the network structure rather 
than a specific network partition thus very convenient to 
estimate the modularity property of the network. 


max tr{Y T G w Y}. (10) 

y t y=i 

According to Rayleigh-Ritz theorem m, the maximum 
for this problem is attained when columns of Y is the 
right eigenvectors U = {u\, ...,ux} corresponding to the 
K largest eigenvalues of the symmetric correlation matrix 
G^. Then the strength of such a community is approxi¬ 
mately equal to its corresponding i-th power of the eigen¬ 
value 


s? 


ulG^u k _ xt uluk _ t 

7 ’ r r / *h 

u k Uk U k Uk 


( 11 ) 


For the convergence of the Potts model across multi¬ 
ple timescales, the vanishing of the smaller eigenvalues as 
the time growing describes the loss of different spin states 
and the removal of the structural features encoded in the 
corresponding weaker eigenvectors. For the purpose of 
community identification, intermediate timescales of local 
uniform states are interesting. If we want to identify z 
communities, we expect to find P t at a given timescale, 
the eigenvalues X\ may be significantly different from zero 
only for the range k = 1, z. This is achieved by deter¬ 
mining t such that |Afc|‘ « 0. 


3. Prove the validity of stability. — Many mea¬ 
sures have been defined to indicate the significance of com¬ 
munity structure, such as modularity Q proposed by New¬ 
man et al |T lTj and spectral cut metrics [Tlj . In [TD] , the 
eigenvector of transition matrix P is also found able to in¬ 
dicate the partition of the nodes in the network. The com¬ 
ponents of eigenvector corresponding to nodes within the 
same community have very similar values and the eigen¬ 
value gaps between different communities can represent 
the significance of the modularity structure. In this part, 
we use phase transition of stochastic dynamic system to 
prove that the stability we proposed is in proportion to the 
eigenvalue gap. Thus, the larger the stability, the larger 
the significance of the community structure. 

Let us demonstrate our argument for the simplest case 
that a network owning 2 communities. As Ai =1, accord¬ 
ing to Eq. we write 


P l — uifi + A 2 U 2 / 2 + | A 3 | f B l , (14) 


where B t is the remainder matrix 


st = E- 

k> 3 




-Ukfk 


(15) 
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Although P need not to be diagonalizable, the repre¬ 
sentation is guaranteed by the nondegeneracy of Ai and 
A 2 . Because Yf,x h( x ) u ^( x ) = 0 and u i > 0, we deduce 
that 

f M = max x<EX f2(x) > 0 > f m = min xeX f2{x) (16) 

Let us now fix a > 0 and according to Eq. m, we consider 
the two nodes sets belonging to two different communities: 

J M («) = {xeX\ fM ~/ M 2{x) < a}, 

I m {a) = {x € X \ < a }, 


The foregoing identity is a fundamental statement about 
the structure of the two communities. From the definition 
of Im(o)iEm{cl) and the identify, we deduce 

E u m ( x ) > 1 - E u m ( x ) < 1 q A2 

x€lIm(o) xGEm(cl) 

(26) 

with analogous inequalities for M replaced by m 
We have seen that 


u 2 {x) 


u m ( X ) - U ln ( X ) 

A !(/m - fm) 


(27) 


These two sets will in fact turn out to be the two phases. 
We also define Em (a) = (pIm{cl ), E m (a ) = ipl m (a) (where 
ip indicates the complement of a set). We take a < 1 and 
note that 


a < 1 =4- Im(ci) and I m (a ) are disjointed. (18) 


Accordingly to the two nodes sets, the corresponding 
the two phases of the system are, 


«mW = M x ) + X t 2 u 1 (x)f M + | A 3 | Bl M , 

i ( x ) = ui( x ) + X t 2 u 1 (x)f m + | A| | Bl m , 


Here B xM (respectively B xm ) is the value of B xy for some 
point y(= y M ) (respectively, y m ) such that f 2 (y ) = /m 
( respectively, f m ). From Ea. (fl9l) . we have 


ui(x) 

u 2 {x) 


) , I \t I /m-B* M - f M 
, fM-fm / I 3 1 , /M-/ m 

Ik . |A,i 

A HfM-fm) + A‘ f M -fm 


( 20 ) 


More generally, for any y one can define 


w* (x) = Ui(x) + X 2 u 2 (x)f 2 (y)+ | A| | B xy = P* y (21) 

We take the scalar product of u l M and with f 2 . There 
are 

Y J fi{x)B t xy = 0 . ( 22 ) 

X 

This follows from 

( 23 ) 

I A 3 I 


and therefore 

^ J f 2 {x)u t M (x) = f 2 (x)u t m (x) = Xy m , (24) 

X X 


From Ea. (l24l) and the fact that )T) m m( ; e ) = X) u m (. x ) = 
1 we can deduce the basic identity 

A 4 ! - A‘ = E «mW(1 - ^r 1 ) = E *&(*)(! - ^) 

x jM x Jm 

(25) 


where r(x) = A| x(linear combination of B xM ,B xm ). 
From the basic identity in Ea. fli^Sl) . we deduce 

Yxx&Im U2 ^ - a[l-(A t 1 -A^)](/M-/ m ) + e ’ / 2 8) 

U ^ X ) ^ a[l-(A*-A*)](/ M -/m) + € ’’ 

where e, e’ are o(l — X\). 

From Eq. KSl) . we notice that the gap of eigen¬ 
vector value is £ x£Im u 2 (x) - £ xeJm u 2 {x) > 

a[1 _ (M _ A ?)]( /M _ /m ) + o(l - A|) and its lower bound 
increases with the stability ©(f) = X\ — X\. Thus, we 
declare that larger stability of communities will extend 
the eigenvector gap between them and thus enhance the 
significance of the community structure. Furthermore, 
one can easily extend the stochastic dynamic system to 
fc-state, k > 2 , in which the stability 0 (f) = Xk — Afc+i 
can also indicate the significance of structure owning k 
communities. 

4. Estimate the spin correlation. — The spin cor¬ 
relation matrix C is very important for the Potts dynamic. 
In this section, we propose a novel way to calculate C us¬ 
ing a new hierarchical block model method based on dif¬ 
ferent granularity (resolution). Stochastic block model [9] 
is a useful tool to detect communities from networks or 
dynamical networks. However, the existing block model 
methods are restricted to the specific task of community 
detection and not suitable to models which need to extract 
multiple levels structure of hierarchical networks. In this 
part, the stochastic block model is extended to a multilevel 
form and exactly coincides with the dynamical process of 
the Potts model. 

Let A nxn be the adjacent matrix of network N , where 
n is the number of nodes. Suppose all nodes of N are di¬ 
vided into L(1 < L < n) blocks, denoted by where 

bu = 1 if node i is in block f, otherwise bu = 0. When 
each block is considered to be inseparable, the granularity 
of network N can be measured by the number of blocks 
g = L. As g decreases from n to 1, the granularity of N de¬ 
generates from the finest to the coarsest. Let B g denotes 
the block matrix B with a granularity g. In particular, 
we have B\ = I nxn . Let matrix Zlxk( 1 < K < L) de¬ 
notes such communities, where K is the community num¬ 
ber and zik = 1 if block l is labeled by community fc, 
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otherwise Zi k = 0. Given Z , define S Kxn, where £kj de¬ 
notes the probability that any node out of community k 
expects to couple with node j; and define C2(uii, ...,u>k ) T , 
where uj k denotes the prior probability that a randomly 
selected node will belong to community k. It is easy to 
show that spin correlation matrix at level g is C g = B g ZS 
so calculating C g is corresponding to estimate B g , Z and 

Define X = ( K 1 Z 1 S, SY) be a pattern unit of network N 
with respect to B g . According to the principle of maxi¬ 
mizing the posterior probability, the optimal X for a given 
network N under B g will be one with the maximal poste¬ 
rior probability. Moreover, we have 


P(X\N, B g ) oc P(N\X, B g )P{X\B g ) (29) 

where P(X\N, B g ), P(N\X , B g ), and P(X\B g ) denote the 
posteriori of X given N and B g , the likelihood of N given 
X and B gi and the priori of X given B g: respectively. 

As discussed above, an optimal X will be the one with 
the maximal P(X\N,B g ) and to maximize P(X\N, B g ) 
is to maximize the product of P(N\X, B g ) and P(X\B g ). 
For a given K, the term P(X\B g ) is a constant, and thus, 
to maximize P(X\N, B g ) is to maximize L(N\X, B g ). 

Let L(N\X, B g ) = lnP(N\X, B g ), and we have 

L K n 

w , b 9 ) = EEEI1 / i a ij ) w k (30) 

1 =1 bi!#0 k—j j— 1 

where f{x, y) = x v {l — x) 1_v . 

Considering the expectation of L(N, Z\X, B g ) on Z, we 
have: 


E[L(N, Z\X, B g )\ 

= Et i Eb ul t o Ek= i 7ifc(E"=i( ln fi&kji a ij)) + Inwfc) 

(31) 

where E[z ik ] = 7 Ik = P(y = k\b = l,X,B g ), i.e., the 
probability of block l will be labeled as community k given 
X and B g . Let J = E[L(N, Z\X, B g )] + \{E k =i w k = 1), 
we have: 


= 0 


dJ 

d J?¥ 

— = 0 

du k 


E[ = lEb .,^0 a ijllk 

W = y?* - — 

2.^1 = i ^o > lk 

Ei=i Eb,,^o 'y ik 

OJk — 

k=i 2^1 = 1 2 ^b U y 6 0 ’ lk 


(32) 


Let P{y = k\v = i) be the probability that node i be¬ 
longs to community k given X and B g , We have: 7 i k = 
P{y = k\b = l,X,B g ) = Eb u ^o T*hb~i P (y = k \ v = 
where is the probability of selecting node i from 

block l. According to the Bayesian theorem, we have: 


P(y = k\v = i) 


P{y = k)P{v = i\y = k) 
E k = 1 p (y = k)P{v = i\y = k) 


(33) 


and 


P(y = k)P{v = i\y = k) = / {£ k j , a,j)w fc (34) 

3 =1 


Thus 


7 Ik = 


n 3 -=l fi£kj> a ij) w k 

Zk 


Si=l ba j .Q Ek= 1 n,=l CLij)w k 


(35) 


As a conclusion, a local optimum of maximizing Eg. (1291) 
will be guaranteed by recursively calculating Eg. (13^1) and 
Eg . (1351) with granularity g. An optimal pattern unit X = 
( K , Z, 5, f2) is calculated given B g and consequently spin 
correlation matrix C g = BgZ’E,. 

For a given network, the hierarchical calculation pro¬ 
cess of C g as g decreases from n to 1 can be incremen¬ 
tally proceeded as follows: First, constructing the ground 
layer by taking each node as one block, and L = g = 
n,B g = J„ X n- Thus, C n = B g ZE = A. Then clustering 
it into n — l communities by selecting a model X n with a 
maximum P(X n \N, B n ). Second, according to X n , form 
B n _ 1 by capsuling each cluster in the ground layer as one 
block. B n _ 1 = B n x Z n _ 1 and C n - 1 = I3 ra _iZ2 = A. 
Then clustering these n — l blocks into n — 2 communi¬ 
ties by calculating a new model X n _i with a maximum 
P(X n -i\N, B n -i). Repeat the second step to construct 
more layers until the process converges, i.e., all blocks are 
grouped into only 1 communities. After calculating all n 
layers C s , the average spin correlation matrix C can be 
taken as the hierarchical average C = (C g ), g = l,..,n. 
One can easily find that the process coincides with the hi¬ 
erarchical dynamical process of Potts model described in 

Figffl 


5. Experiments. — To show that the model can 
uncover hierarchical structures in different scales, FigE] 
and Fig J3] give two examples of the multi-level community 
structures, RB 125 network m and H 13-4 network [3]. In 
both examples, the most persistent A reveals the actual 
number of hierarchical levels hidden in a network. The 
significance of such levels can be quantified by their corre¬ 
sponding length of persistent time. Longer the time per¬ 
sists, more robust the configuration is. From Fig |2(b)| and 
Fig 3(b)] we can observe 25 and 16 are the optimal num¬ 
bers of communities in RB 125 and H 13-4 networks owning 
the longest persistence, respectively. However, 5 modules 
and 4 modules are also reasonable partitions which show 
another fuzzy level of the hierarchical networks. These 
results are perfectly consistent with the generating mech¬ 
anisms and hierarchical patterns of these two networks. 

Furthermore, we also show that the variation tendency 
of stability 0 (r) in the two cases shed a light on the spin 
configuration. From Fig ]2(b)| and Fig ]3(b)[ there are some 
local maximal values representing better community struc¬ 
ture. Thus, we can find these local maximal timescales T 
corresponding to the desirous number of communities and 
apply G T to a specific partition method. Furthermore, the 
stability will reach the lowest value at the end time of all 
A. The stability begins to increase when it transits to a 
new state. One can use O(r) to estimate the modularity 
property of complex networks, and the larger the 0 is, the 
stronger the network community structure. So, one can 
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Fig. 2: (a) Structure of RB 125, with 25 dense communities and 
5 sparse communities, are highlighted in the original network, 
(b) The value of A(r) and @(t) versus time t. 


Fig. 3: (a) Structure of H 13-4, with 16 dense communities and 
4 sparse communities, are highlighted in the original network, 
(b) The value of A(r) and 0(r) versus time t. 


find the largest corresponding 0 value for a specific num¬ 
ber of community A and use it to indicate the robustness 
of modularity structure. For H 13-4 shown in Fig |3(b)| 
the stability of 16 communities structure, r( 16 ) = 0.62 
when t = 3, is larger than T(4) = 0.43 when r = 12. 
This indicates that the community structure containing 
16 modules is more robust than community structure con¬ 
taining 4 modules. Similarly, for RB 125 network shown 
in Fig |2(b)| T(25) = 0.71 corresponding to 25 communi¬ 
ties structure when r = 3 is larger than T(5) = 0.18 when 
r = 13. The robustness of community structure indicated 
by soft stability 0 favors finer but obvious modules which 
reasonable for many real networks. In addition, the dif¬ 
ference between the stability measure we proposed and 
the modularity Q 2] is emphasized. We also applied our 
framework to the hierarchical network with different mod¬ 
ular sizes and some representative real networks. Finally, 
the relationships between our work and some famous con¬ 
cepts proposed in [6] and 0 are analyzed. These results 
are shown in the part of Supplementary Material [I?. . 

6. Conclusion. — In summary, we have presented 
a more theoretically-based community detection frame¬ 
work which is able to uncover the connection between 
network’s community structures and spectrum properties 
of Potts model’s local uniform state. Important informa¬ 
tion related to community structures can be mined from a 
network’s spectral significance through a Markov process 
computation, such as the stability of modularity structures 
and the optimal number of communities. Our method 
does not provide a unique optimal partition for the graph. 
Rather, we obtain number of stable levels and stability at 
each level over different layers of the hierarchical struc¬ 
ture. Its effectiveness and efficiency have been demon¬ 


strated and verified both theoretically and experimentally. 
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